Anomalous Light Absorption by Small Particles 
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A new type of resonant light absorption by a small particle (nanocluster) is reported. The problem 
cannot be described within the commonly used dipole scattering approximation and should be 
studied with methods based upon the exact Mie solution. It is shown that the absorption cross- 
section has giant maxima realized at small values of the imaginary part of the complex dielectric 
permittivity of the particle. The maxima are situated in the vicinity of the plasmon (polariton) 
resonances and correspond to the regions where the dissipative damping equals the radiative one. 
The case is similar to the recently introduced anomalous scattering [PRL 97, 263902 (2006)] and 
exhibits similar peculiarities. 

PACS numbers: 42.25.Bs, 42.25.Fx, 78.67.Bf, 78.20.Bh 



Since the first quantitative study by Lord Rayleigh in 
1871 [l[, the problem of light scattering by small particles 
has remained one of the most important and appealing 
issues of electrodynamics. There are thousands of articles 
and numerous monographs devoted to this matter, see 
e.g. [H and references therein. Nonetheless, the problem 
is far from completion. 

Recently it was pointed out [1, 0, [sf that small parti- 
cles made of weakly dissipating materials may exhibit the 
so-called anomalous scattering (AS), which has very lit- 
tle in common with Rayleigh scattering, other than both 
occur without changes in the frequency of incident light 
00. AS results in giant optical resonances with narrow 
linewidths and inverted hierarchy, i.e., the partial extinc- 
tion cross-section for the quadrupole resonance is larger 
than that for the dipole one, etc. [1, 0- It occurs when 
the dissipative losses are small enough to become negligi- 
ble relative to the radiative damping [1, 0, 0]. Therefore, 
it seems that the phenomenon has nothing to do with the 
maximal absorption of the light by the particle, when the 
dissipation should be as large as possible. However, this 
is not the case. 

In this Letter, the absorption cross-section 
(Jahs{R-,^'p-,^p) of a spherical particle with radius R 
and complex dielectric permittivity = + ie'^^ is ana- 
lyzed, based upon the exact Mie solution [8]. The exact 
solution automatically includes both types of damping 
(dissipative and radiative). Thus, the problem is just to 
find an appropriate asymptotic of the cumbersome Mie 
solution. Such an analysis indicates that (7^55 (i?, e^, e^') 
for small particles exhibits sharp giant peaks whose 
maxima lie at e'^ <c; 1, while the corresponding values 
of e'p coincide with those for plasmon (polariton) res- 
onances. It should be stressed that apart from purely 
academic interest, the problem of maximization of cFahs is 
important for a number of applications of nanoparticles 
in technology, biology, medicine, etc. - which provides 
additional stimuli for the present study. 



To understand the physical nature of the aforemen- 
tioned peculiarities of (Jahsi let us estimate the power P 
dissipated in a unit volume of the particle. Neglecting 
the magnetic part of the dissipation owing to the parti- 
cle smallness [9[, one obtains P ex e''|Ep, where E(r) is a 
complex amplitude of the monochromatic electric field in 
a given point inside the particle. On the other hand, in 
the vicinity of each plasmon (polariton) resonance E(r) 
undergoes a sharp increase. At these resonances in the 
case of Rayleigh scattering, E(r) is cut off by dissipative 
losses, so that E(r)£ oc Thus at the resonances, 

P (X 1/e^' ^ 00 at e^' ^ 0, which gives rise to the corre- 
sponding divergence of cFahs- 

The divergence occurs because the approximation used 
is invalid. In fact, (Tabs does not diverge. At small enough 
e^', the radiative damping prevails over the dissipative 
losses and Rayleigh scattering is replaced by AS, which 
provides a different cutoff for E(r) [3, 0, S 0]- 

To inspect the problem accurately, one should em- 
ploy the exact Mie solution. It is convenient to 
introduce dimensionless cross-sections Qext,sca,abs = 
(7ext,sca,ahs/^R^, whcrc CFext, CFahs and Gsca respectively 
stand for the extinction, scattering and absorption cross- 
sections of the particle. Then, the Mie solution yields the 
following well-known expressions: 
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The summation is over the corresponding partial multi- 
pole cross-sections, the so-called size parameter q equals 
TimRko^ rim = V^rn Purely real refractive index of the 
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environmental medium, and ko stands for the wavenum- 
ber of the incident hght in vacuum {ko = uj/c^ where c is 
the speed of hght). 

Amplitudes ai and bi may be presented as follows [3[: 
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Here e = ep{u;) / em{^) and 
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Cf^ for terms G\^^ and —F^^^ in the denominator of Eq. (pT]) 



' ^ are expressed in terms of the Bessel 
[Ji+i/2{^)] and Neumann [A/'^+i/2(^)] functions, whose 
expansions in power series gives rise to the following for- 
mulae, valid at small q\ 
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(ellipses denote omitted higher order terms in q). The 
explicit expressions for F^^^ and G^^^ are not required 
because of the estimate \bi\ ^ q^^^^ <C |a^|, which allows 
neglect of hi relative to at g <C 1. 

Eq. (|4|) results in the following expression for the par- 



tial absorption cross-sections Q^^^g ~ 
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where prime and double prime denote the corresponding 
real and imaginary parts, respectively. 

The plasmon (polariton) resonances are defined by the 
condition Gi^\{q, e') = 0, which determines the resonant 
values of the real part of the dielectric permittivity e[{q) 
and through the dependence e(cc;), the resonant frequen- 
cies ooi{q). For the problem in question the maxima of 
aabs should be situated in the vicinity of the resonances. 
Let us inspect the vicinity of a certain ith resonance, 
presenting the dielectric permittivity in the form 



e = 
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where 5e' and e" are small, see Eq. (|6j). Then, the leading 



approximations for G^f^ and Q^^^ read: 
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According to Eq. (j9|), F^^ is negative and does not 
depend on (5e', ^' . It allows one to transfer from (5e', ^' 

to G^f^ , G^f^ , considering the latter as new indepen- 
dent variables, see Eq. (p!Q|) . The maximum of Q^^^ cor- 



responds to G 



(i.e., it is achieved exactly at the 



(4) resonant frequency uji) and G 



[ay 



-F^""^' [it should be 
remembered that G\^^ > 0, see Eq. (fTQ|) ]. Note, the two 
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correspond to the dissipative and radiative damping, re- 
spectively. So, is maximized when the former equals 
the latter. 

Thus in the vicinity of each plasmon (polariton) reso- 
nance, Q^^^^ passes through the local maximum 

2^ + 1 
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which equals 1/4 of the maximal AS resonant partial 

extinction cross-section, Q^^xtmax 01' ^^^^ linearly 
increases with an increase in the order of resonance I. 
However, in contrast to Q^f^t^ax (achieved at e'' = 0), 



the maximum of Q^^^ is realized at 
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inewidths n 
J 2 are as follows: 
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The corresponding linewidths measured at the level 
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Though the values of Ae', Ae'' are close to each other, 
the lineshapes are quite different. To study this issue in 
detail let us introduce new variables: 
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Then, (in the given approximation) the resonant partial 
absorption cross-section is reduced to the following uni- 
versal, q- and ^-independent form: 
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(18) 



Thus, the line along ^ (i.e., e') axis has a typical sym- 
metric Lorentzian shape, while the one along ( (i.e., e") 
axis is strongly asymmetric, see also Fig. [H 

It is interesting to mention that at the resonant points 
of the absorption cross-sections given by Eqs. (p!3|) - (p!4|) . 
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FIG. 1: (color online). Resonant absorption lineshapes along 
axes of real (a) and imaginary (b) parts of dielectric per- 
mittivity of the particle, in rescaled dimensionless variables 
Eq. ()17p . The universal profiles corresponding to the approx- 
imate Eq. ([18)) are shown as dotted lines. Lines for the dipole 
and quadrupole resonances derived from the exact Mie solu- 
tion for q — 0.3 are shown as solid red and blue lines, re- 
spectively. Note high accuracy of the approximate expression 
Eq. ([18]); while for the dipole resonance the universal lines 
given by Eq. ()18p are still a bit narrower than those for the 
exact Mie solution, for the quadrupole resonance they are 
practically identical with the latter. 



the scattering cross-section is Qica(4'40 



^ abs max'' 

see Eq. ([T2]) . At this point the partial extinction cross- 
section equals just double of Q 



sea max • 



The dimensional absorption cross-section which cor- 
responds to Eq. m <^al ma. = 7r(2^ + l)/2(n„fco)2 
does not vanish at ^ 0, which looks confusing. The 
confusion is resolved if one notes that for the function 
^a6s(^' e^', e/')^ the point R = is singular; the function 
does not have any limit at this point and may take any 
value varying from to a'^^^ depending on the way 
to approach this point in (i?, e/ ^ e^") space. In particu- 
lar, if one first fixes values of e', e" and then tends R to 
zero (which is the closest case to a possible experimental 
situation), cr^^^^[R^e\e") vanishes, as it should be for a 
particle with zero radius, cf. the analogous properties 
of the resonant partial cross-sections of the anomalous 
scattering 0,0,11. 

Regarding off-resonant partial cross-sections, their val- 
ues are given by the expression: 
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Comparison with Eq. (fT2]) reveals that at small q in the 
vicinity of the resonances, the contribution of the par- 
tial resonant cross-section to the net cross-section Qabs 
is overwhelming. 



Let us discuss possible experimental observation of the 
anomalous absorption and related issues. In the pre- 
ceding theoretical analysis three quantities i?, e' and e" 
have been regarded as independent. In reality it is dif- 
ficult (if possible) to tune e' and e" independently. Ac- 
tual independent control parameters are R and uo. Then 
to maintain the resonance conditions, these quantities 
should satisfy the following set of equations: 

e'{uj,R) = e[{uj,R), 
e"{uj,R) = e'l{uj,R), 

where the left-hand sides follows from the dispersion 
properties of the particle material (including spatial dis- 
persion) and the right-hand sides are given by Eqs. ([T3|) . 
([T4|) . respectively. Solutions of the equations (if any) 
yield a unique discrete set of pairs {uoi^Ri)^ where Ri 
should satisfy the additional constraint following from 
the restriction q <^ 1. These conditions are very strict 
and the entire set of them is very difficult to fulfill in 
any real experiment. Thus, it seems the phenomenon 
discussed is just a "virtual" effect. 

Fortunately the situation is not so dramatic. There 
are several reasons for that. First, the resonance lines 
are narrow in the (e', e") plane, but this is not necessarily 
the case for the linewidths along the uo axis. The latter 
is determined by the dispersion properties of the particle 
materials. For most materials in the discussed region 
(e' < —1), the dispersion of e is rather weak and the 
resonance conditions imposed on uo are not as strict as 
those for e. Second, while the decay of the absorption 
cross-section at e' < is very sharp, it is rather slow 
at > e^, see Fig. {[IJd). Finally, one has to take into 
account sharpness of the i?-dependence of the right-hand 
side of the resonance condition, see Eq. (fT4|) . Due to 
all these facts, to observe the anomalous absorption it 
suffices to select material with small values of e" for e' 
lying in the region of possible resonances (in reality for 
—2.5 < < —1.5). Then, varying R (or/and uo), one 
inevitably passes through the vicinity of the resonance (s). 

To illustrate this general reasoning, the absorption 
cross-section of aluminium nanoparticles is studied in de- 
tail. The model employed in the study is identical to the 
one discussed in Ref. [lO[: first, the dependence of e on 
uo is taken from an empirical table [iT]. Then, it is ap- 
proximated by the Drude formula: 
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(20) 



To enhance the accuracy of the approximation, quantities 
uop and 7 are regarded as functions of uo, calculated at 
every point in the table. This allows the identity of the 
initial table with the one resulting from Eq. ([2Q|) . To take 
into account collisions of free electrons with the particle 
surface, the obtained 7(0;) for bulk aluminium is replaced 
by Jeff = 7(^) -\-vf/R, where vp = 10^ cm/s stands for 



4 



the Fermi velocity of the free electrons. Next, polynomial 
interpolation between the table points is used to get a 
smooth e(a;), which is used to calculate the absorption 
cross-section derived from the exact Mie solution. 
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FIG. 2: (color online). Contour plots of the dimensionless 
(a) and dimensional (b) absorption cross-sections for an alu- 
minium particle calculated based upon the exact Mie solution 
and the actual empirical dependence e(a;) for aluminium [ll[. 
Note that the localized maxima corresponding to different 
orders of the resonances for the cross-section normalized over 
TiR^ (a) become more extended and merge for the dimensional 
cross-section (b). 



TABLE I: Approximate theory vs. exact Mie solution for Al 



i 


A^, nm 


Ri, nm 
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T 


Q 


Appr. Exact 


Appr. Exact 


1 


143 


6.45 


0.283 


-2.16 


-2.19 


0.046 


0.216 


0.617 


2 


136 


19.27 


0.890 


-1.67 


-1.88 


0.047 


0.159 


0.901 



Results of these calculations are presented in Fig. [2l 
Quantitative comparison of the results with the devel- 
oped approximate analytical theory is shown in Table H 
The resonant values Xi and Ri correspond to the local 
maxima of Fig.[T^. The quantity is calculated as a root 
of equation G^'^^^(<7, e') = 0, see Eq. (|6|). The quality fac- 
tor T equals the ratio of the net absorption cross-section 
(at a given for the aluminium particle) to Q^^^ 
see Eq. (fT2|) . 

Note that the obtained e'^ for the aluminium particle 
are considerably larger than those which, according to 
Eq. ([Hj), should correspond to the resonances; neverthe- 
less well-pronounced maxima of the dipole, quadrupole 
and even octupole resonances are seen on Fig. [T] straight- 



forwardly. The absolute maximum of the relative ab- 
sorption cross-section, MaxlQabs} ^ 11.516, is achieved 
at A 143 nm; R ^ 6.45 nm and corresponds to the 
dipole resonance {£ = 1). 

The presented analysis has thusly revealed a small par- 
ticle made of weakly dissipating materials has the absorp- 
tion cross-section which may exceed its geometric cross- 
section in order of magnitude, and even more. Though 
the anomalous absorption and the anomalous scattering 
correspond to opposite limits (the maximal and vanishing 
absorption, respectively), they have much in common, 
namely both are associated with giant amplification of 
incident light inside the scattering particle in the vicin- 
ity of the plasmon (polariton) resonances accompanied 
by interplay of the radiative and dissipative dampings, 
both cannot be described within the framework of the 
dipole approximation and require analysis of the exact 
Mie solution, and both may exhibit similar unusual prop- 
erties (the inverted hierarchy of the resonances, singular 
dependencies of cross-sections on the particle size, etc.). 
The calculations for aluminium particles indicate that 
the anomalous absorption should be an experimentally 
observable phenomenon. A future endeavor will be to 
carry out the corresponding experimental study, which 
hopefully may be motivated by this Letter. 
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